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We provide software with a graphical user interface to calculate the phenomenology of a wide class 
of dark energy models featuring multiple scalar fields and potentials that are arbitrary functions 
of exponentials. The user chooses a subclass of models and, if desired, initial conditions, or else a 
range of initial parameters for Monte Carlo. The code calculates the energy density of components 
in the universe, the equation of state of dark energy, and the linear growth of density perturbations, 
all as a function of redshift and scale factor. The output also includes an approximate conversion 
into the average equation of state, as well as the common (wo,w a ) parametrization. The code is 
available here: http: //github. com/kahinton/Dark-Energy-UI-and-MC 


I. INTRODUCTION 

The puzzle of dark energy is one of the most important 
outstanding questions in all of physics. Studying the ob¬ 
servational signatures of dark energy - and using those 
signatures to understand the nature of dark energy - has 
been a highly active area of research since the first signs 
of dark energy were discovered over 15 years ago. (For a 
review, see e.g. Frieman et al. [d].) Dark energy models 
have a rich phenomenology, leaving varied fingerprints 
on the growth of cosmic structure and on geometrical 
quantities in the universe. 

We have written a piece of software that makes it 
easier to compute quantitative predictions for dark en¬ 
ergy models broader than a simple cosmological constant. 
These models possess a richer phenomenology that can 
be tested with ongoing and future experiments. They 
also complement the choice of dark energy models - and 
computer code - made by others [2]- Our software is 
capable of handling a wide class of models with multi¬ 
ple scalar fields, and specifically makes two subclasses 
particularly easy to calculate with a simple user input. 
In these models, the user can calculate the expansion 
history and the growth of structure in the linear regime, 
therefore enabling the user to obtain various implications 
for observable quantities in cosmology. The code also re¬ 
turns an approximate conversion of the equation of state 
of dark energy w(a) into two commonly used parameters 
(wo,w a ), as well as several other commonly used param¬ 
eters. 


II. DARK ENERGY MODELS 

We choose two classes of “assisted dark energy” mod¬ 
els, based specifically on the analysis by Ohashi and Tsu- 
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jikawa [3], and originally developed by Refs. HE]. The 
starting assumption is the existence of n scalar fields, and 
a Lagrangian density that takes the following form 

n 

Aelds = J2x ig (X ie x ^) (1) 

i= 1 

where Xi = is the kinetic term for i-th 

field, and g(Y) is an arbitrary function of its argument 
Y = X e A( ^, where A is a dimensionless parameter. The 
existence of scaling solutions dictates that this be the 
most general form of the Lagrangian density for expo¬ 
nential potentials Hog. a noteworthy special case of our 
class of models is quintessence with a single field (n = 1), 
and a simple exponential potential V{(j)) m cexp(A0), for 
which g(Y) = 1 — c/Y. As we explain below, this class 
of models also includes ghost condensate models [6] fea¬ 
turing a background scalar field which, despite having a 
kinetic term of the “wrong” sign, are stable and feasible 
due to suitable higher-order kinetic terms. 

Our code allows a user to test models with an arbitrary 



FIG. 1. Sample screenshot of the software, showing text in¬ 
puts for dark energy model parameters, as well as checkboxes 
for selecting the desired plots. 
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FIG. 2. Example plots generated by our program, showing output for a single model for our quintessence (left panel) or ghost 
condensate (right panel) class of models. The user can choose to plot this for any arbitrary choice of model parameters, as a 
function of either scale factor a or redshift z. 


function g(Y ), for an arbitrary number of fields n. As an 
example we also provide two special cases which have 
been run through the Monte Carlo generator, and whose 
results can be accessed with very simple inputs in the 
graphical user interface (GUI): 

g(Yi) = 1 — Ci/Yi (quintessence) (2) 

g(Yi) = — 1 + CjYi (ghost condensate). (3) 


In each case, the temporal evolution of all quantities is 
most easily tracked in terms of scaled variables 
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where p r is the radiation energy density. The evolution 
equations, expressed in terms of the time variable N = 
In (a), are [3 
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We numerically evolve these equations. Then all phys¬ 
ical quantities of interest can be found, for example the 
energy densities of matter and dark energy 


Cde = ^iXi) + 2y iff 'c m (7) 

i =1 


and the dark energy equation of state 
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where the time-dependent parameter Y{ can be calcu¬ 
lated via Yi = x^/g?. Note that for quintessence, the 
functional form of g(Yi) in Eq. © implies that the po¬ 
tential is the sum of exponential potentials for each field, 

V{<t>) = J2i °i exp(Aj^j). 


III. USING THE PROVIDED SOFTWARE 

Our code can be found at the following 
GitHub repository: http: //github. com/kahinton/ 
Dark-Energy-UI-and-MC. The code gives the user 
several tools for testing and analyzing various models 
of dark energy. The first step in being able to use 
these tools is to run a model through our Monte Carlo 
generator. The user must supply several inputs to the 
MC generator: a specific function g(Y), as described 
above; a name to differentiate the model from others; 
and the number of different initial conditions to test the 
model over. The provided name will be used to label 
this model in the user interface. The user must also 
specify the values for several parameters governing each 
model. Specifically these parameters are the number 
of fields that will be acting as dark energy, n, as well 
as the range for the initial conditions of yi, q, and 
u. In our code we adopt a simplification suggested in 
[3], which assumes that a single field will act to assist 
inflation in the early universe before becoming small, 
while the remaining fields will be small in the early 
universe before becoming the dominant form of energy 
at late times. The “assisted inflation” field will have 
one set of unique initial conditions, while the remaining 
n — 1 fields will all share the same initial conditions. 
This later simplification drives the value of wde as close 
to —1 as possible at late times, and also allows models to 
be tested at the same speed independent of the number 
of fields being used. 

By using these simplifications the user only needs to 
specify the range of initial conditions for the first two 
fields. The code then assumes, for example, that all val¬ 
ues of X 2 to x n are equal. To ensure that a scalar field 
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FIG. 3. Left panel: the first two principal components used to convert from w(z) to w ave and (wo,w a )', note that w ave requires 
only the first principal component. The components are based on the NASA/DOE Figure of Merit Science Working Group 
analysis 7j; see Appendix for details. Right panel: a sample scatterplot made from the user interface. Here one can see the 
results of wq and w a for 10,000 pre-supplied ghost condensate models, case with 15 fields. 


dominated solution is produced, the values for A^ are 
chosen automatically, according to the inequalities (see 
Eqs. (29) and (32) in [3|) 

— > 88.9; Am < 2 , (9) 

P,x P,x 

Here p is the Lagrangian density, and p,x = dp/dX. 
When the code is run, for each specific test, the initial 
conditions for each parameter are chosen from a random 
flat distribution over the selected range. When all of 
the tests have finished running the model will be added 
to the user interface automatically. As an example of 
what is output to the user interface, we have included 
the results of tests of quintessence and ghost condensate 
models having 5, 10, 15, and 20 fields. 

The user interface also contains two optional tools. 
The first of these allows the user to run a single instance 
of any model having been run through the Monte Carlo 
generator. Each of the required parameters are assigned 
a set of fiducial values that depend on which model is cho¬ 
sen; however, the range over which the model has been 
tested is also supplied, allowing the user to experiment 
with values other than those provided. Once the user se¬ 
lects values for these parameters, they can choose to plot 
various quantities as a function of either scale factor a or 
redshift z. These quantities include: 

• Energy densities in units of critical (radiation 
D#(a), matter f1 m(&), dark energy Hde(&)); 

• Equation of state of dark energy w(a ); 

• Linear growth of density fluctuations D(a ) (we plot 
the quantity D(a)/a). 

In order to provide the growth of linear perturbations, 
concurrently with the other equations we also evolve 

D" + D\H'/H+2)-3/2(Ho/H) 2 n M (l+z) 3 D = 0 (10) 


where D(a ) = (Sp/p)/(Sp/p ) a= i is the linear growth rate 
of matter perturbations, and where the primes indicate 
a derivative with respect to ln(a). The quantity D/a 
for the selected model will immediately be displayed on 
screen for the user to examine. 

The second option available to the user is to examine 
the results from the Monte Carlo generator. Specifically, 
this tool allows the user to generate two-dimensional scat¬ 
terplot s on the fly for further analysis - for example, 
to see what range of phenomenological outcomes is pro¬ 
duced by a given model or class or models. Each point in 
a given scatterplot corresponds to an output of a single 
model. The output can be any of the following parame¬ 
ters or functions: 

• Initial conditions from each model including aq, x 2 , 
2/i, 2 / 2 , ci, c 2 , Ai, A 2 , u , and n; 

• Energy densities in matter and dark energy today 
VLm and Dde, as well as the equation of state today 
WBE{a = 1 ). 

• Effective “average” equation of state re ave , calcu¬ 
lated from the first principal component, as de¬ 
scribed in the Appendix. 

• Effective values for the parameters wq and w a based 
on parametrization [8] w(a) = wq + w a (l — a), de¬ 
rived from the first two principal components as 
described in the Appendix. 

An example of the Monte Carlo output showing the wq — 
w a plane is shown in the right panel of Fig. [3] 
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Appendix: Principal components, w ave , and ( wo,w a ) 

Here we explain in more detail how to calculate the 
aforementioned quantities wo, w a , and w p from the prin¬ 
cipal components of the equation of state, following 
mm- First, we formally expand the equation of state 
in terms of principal components 

1 + w(a) = X a»e*(a), (11) 

i 

where aq is the coefficient of the i-th principal component 
ei(a). The idea is to convert the first principal compo¬ 
nent cti into the averaged value of the equation of state 
re ave , and first two principal components (oq,^) into 
(wo,w a ). This approach is justified because the first few 
principal components carry essentially all the necessary 
information about the effects of dark energy dynamics on 
the expansion of the universe on observable scales. 

The principal components e$ can be determined for any 
given dataset; here we use the Figure of Merit Science 
Working Group’s (FoMSWG) publicly available code [7], 
and a combination of Planck+BAO+SN+WL available 
in their code. Note that, while the component shapes 
— especially the all-important peak value of e± (a) that 
shows the temporal epoch of maximum sensitivity to the 
equation of state — depend on the cosmological probe as 
well as specifications of a given experiment, once we com¬ 
bine the probes the pull of different probes is expected 
to average out, leading to a fixed set of shapes. 

While the normalization of the e*(a) is arbitrary in 
principle, the FoMSWG principal components that we 
use are normalized as f ef(a)da = 1. The coefficients ol{ 
can be obtained as 

ai = J [1 + w(a)] ei(a) da ( 12 ) 

where w(a) is the actual equation of state of the theoret¬ 
ical dark energy model we are studying. 


The final step is converting the first principal compo¬ 
nent into re ave , and the first two into wo and w a ; we do 
this via [10| 


1 + Wave = -y (average w) (13) 
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where cti are defined in Eq. (12), and 



(14) 


(15) 


Equations (13) and ( |14| ) are now our definitions of the 
parameters w ave and (rco,w a ), respectively, given a dark 
energy history w(z) which determines the aq. Previous 
work m confirms that the two-parameter equation of 
state closely follows the true w{z) over the redshift range 
most effectively probed by the data. 

Note too that the constraint wq > —1, which follows 
from w(z) > —1, is not strictly obeyed by wo obtained in 
this way since wo and w a are now essentially a fit to the 
dark energy equation of state history. In other words, 
the derived parameter wq can be slightly smaller than 
— 1; this is again because we are fitting these three re- 
parameters and forcing them to the model’s actual equa¬ 
tion of state evolution. These parameters are neverthe¬ 
less provided since 1) they provide a very useful test bed 
to gauge the effectiveness of the range of dark energy 
model behavior, and 2) they are fairly representative of 
these models, since actual dark energy models often do 
exhibit the wq + w a (l — a) scaling. 
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